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Abstract. Polymers in confined spaces lose conformational entropy. This induces 
a net repulsive entropic force on the walls of the confining space. A model for this 
phenomenon is a lattice walk between confining walls, and in this paper a model of an 
adsorbing partially directed walk is used. The walk is placed in a half square lattice 
with boundary 9L^, and confined between two vertical parallel walls, which are vertical 
lines in the lattice, a distance w apart. The free energy of the walk is determined, as 
a function of w, for walks with endpoints in the confining walls and adsorbing in dh\. 
This gives the entropic force on the confining walls as a function of w. It is shown 
that there are zero force points in this model and the locations of these points are 
determined, in some cases exactly, and in other cases asymptotically. 
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A linear polymer in a good solvent placed near a geometric obstacle (such as a hard 
wall) loses conformational entropy. This loss of entropy induces a net force on the 
obstacle 1[^|, and if the polymer can move freely, then it will tend to move away from 
the obstacle. The induced force is repulsive, and the polymer exerts an average net 
pressure on the wall. These entropic pressures have been seen in experiments 
and have been modelled numerically (see for example reference [^). 

A similar situation is seen when a polymer is placed in a confined space; the loss in 
conformational entropy induces a net repulsive force on the walls of the confined space. 
For example, a polymer between two colloidal particles loses entropy as the particles 
approach one another; this induces a net repulsive force between the two particles. The 
repulsion between particles is the mechanism underlying the stabilization of a colloid by 
a polymer |[^|. This phenomenon was examined numerically using a self-avoiding walk 
model of a polymer confined between two hard walls |[^|. A self-interaction self-avoiding 


walk model of polymer pulled at its endpoints was similarly considered in reference |13,21 
See references l27l[30l[3T1 for more results. 


In this paper a two dimensional partially directed walk model is used to model the 
entropic forces induced by an adsorbing polymer between confining walls. In addition, 
the entropic pressure of an adsorbing polymer on the adsorbing wall will be examined 
as well. Similar models were examined in reference 1[T7]|, using directed path models to 


model the pressure of a directed path on the adsorbing boundary. Related work on the 
entropic pressure near knotted lattice ring polymers was reported in reference |[^|, which 
built on the results in references 


In this paper we continue the study of entropic pressure near lattice polymers by 
using a partially directed walk model in the square lattice. Directed path models were 
examined in reference |[T^, where the entropic pressure near directed paths and staircase 
polygons models (of a grafted linear polymer) was calculated. In particular, for a directed 
path of length n from the origin in the positive half-lattice, the pressure on the x-axis a 
distance proportional to an from the origin is given by 

2V2 


P.(a) = 




TT n^ a^ 


(1-a) 


+ 0(n-^/2) 


This, in particular, shows that 
lim n^/^P,{a) = 

n—^oo 


2 V 2 


V^TT a3(l-a) 

as n — 00 . 

Partially directed walk models of polymers were introduced in references 
have been widely used as models of polymer entropy; see for example references |4 


( 1 ) 


( 2 ) 


and 
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The generating function of adsorbing partially directed walks has been computed in 
several models using the Temperley method j^|, and in this paper we will follow a 
similar approach to that of reference [^ |. 

Let be the square grid with standard basis {ei, 62 }. A point or vertex x G Z^ 
is a vector with Cartesian coordinates (x(l),x(2)). Two points x,y G Z^ are adjacent 
in Z^ if ||x — y ||2 = 1. The edge (x~y) is a unit length line segment between adjacent 
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Figure 1: A partially directed walk from the origin in of width w and endpoint at height h. The walk is 
weighted by the activity a conjugate to the number of edge-visits in 9L^ and its endpoints are tethered to 
two vertical walls a distance w apart which bounds a confining space (vertical slit of width w) containing 
the walk. If h = 0 then this is an adsorbing bargraph path from the origin in L^. 


points (or vertices) x and y in I?. The square lattice is the set of all edges between 
adjacent vertices of defined by 

= {(x~y) I x,y eIF and ||x—y ||2 = 1}. (3) 

The positive square lattice is given by 

= {(x~y) G I x(2) > 0 and y(2) > 0}. (4) 

The boundary of is given by 

= {(x~y) G I x(2) = 0 and y(2) = 0} (5) 

and every edge in is parallel to ei and normal to 62 . 

A directed edge is a directed unit length line segment from a vertex x to a vertex 
y in It is denoted by {x^y). A directed edge {x^y) is in the East direction 
if x(l) + l = y(l) and x(2) = y(2), it is in the South direction if x(l) = y(l) and 
x(2) —1 = j7(2) and it is in the North direction if x(l) = j7(l) and x(2) + l =y(2). 

A partially directed walk cu in of length n is a sequence of n directed edges in the 
East (E), North (N) and South (S) directions, starting in 0, such that a directed edge 
in the N direction cannot immediately be followed by a directed edge in the S direction 
and vice versa. That is, cu is a sequence of directed edges in the E, N and S directions 
((v^^vj), (vj'^'i^). {Fz^v-Fj ,..., {vF-i^Vn)) such that Vq = 0 and all the Vj are distinct 
(so that u is self-avoiding). The terminal or last vertex of u is Vn. 

If a directed walk from 0 is in the positive square lattice L^, then it is a positive 
partially directed walk. A positive partially directed walk is illustrated in figure [1} 

A partially directed walk from the origin in the positive half-lattice is illustrated in 
figure[^ Edges in the walk which are in the x-axis are weighted by the generating variable 
a. For large a the walk stays near the x-axis, so that this is a model of an adsorbing 
partially directed walk. In section |^the generating function of walks of fixed width w 
is determined using the kernel method (^j^. In particular, we show that the generating 
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function of partially directed walks from the origin in the half-lattice, of width w with 
final vertex at height h, is 


W W 

s=0 k=0 ;>0 


W+l — S 
W+i+l — S 


w — k — s 




y/7+2/(^_l)/c 


( 6 ) 


The variable y generates vertical steps (and there are exactly w horizontal steps). This 
result can be used to extract the partition function of walks of length n, width w and 
last vertex at height h (see equation (37)). 

In sectionj^the entropic forces in the model are examined in the scaling limit, defined 
by putting w = [ctnj for walks of length n, and then taking n —oo in the model (see 
figure [^. This gives the limiting free energy F{oc) of the walk, which is computed for 
several cases. For example, for walks with h = 0 and a = 1, the limiting free energy is 
given by 


Fo(ct) = a log(2cf(l —Cf)) — {Vc+2a — l) log('\/C+2ct —1) 

- [1-VC) log(l--/C) + ^VC log (^ 


C±a\ 

C-aJ' 


(7) 


where C = 2a'^ — 2a + 1. The derivative of this expression to ol gives the limiting force 


•^o(a) = ^Fo{a) = log 


/ 2a(l-a) 
UC + 2a-l)2 


+ ^ log 


(C-Q:)(2a+C-1)2 

(C+c.)(l-C)2 


( 8 ) 


It follows directly from the above that T'o(ct) = ——ct) (so that To(|) = 0). That is, 
the limiting force vanishes (in the limit) when the confining walls are a distance w = 
apart. This is the zero force point in the model, and it is located in this model when 
the walk is extended one-half of its length along the adsorbing line (or when one-half 
of the edges (steps) in the walk is in the horizontal direction). The limiting force curve 
JFo(q:) is also symmetric on reflection through the point (|,0) in the (ct, T'o(Q:))-plane. 
For small ct > 0 the above shows that T'o(ct) ~ | logcf| (see equation (|6^). 

In addition to the above, the asymptotic forces in other models are examined. These 
models include walks with a lifted endpoint, as well as walks which are adsorbing. In each 
case asymptotic expressions for the limiting forces are determined, and the location of 
the zero force point is determined. 

In section Infinite size asymptotics of the model is developed for the case that a = 1. 
The partition function simplifies to a single summation which may be approximated using 
a saddle-point method. We determine an asymptotic expression for the partition function, 
and we use this expression to determine the asymptotics for the repulsive force between 
vertical walls (see section [^. In particular, if the walk has length n and the confining 
walls are a distance an apart, then the asymptotic force has leading term asymptotics 
fFn{a) = —^ + (n — |) — n log ct + 0{na) as ct \ 0. For small values of a the force is 
positive (repulsive). 

In sectionj^the model is examined numerically. This is done in particular to (1) verify 
the asymptotic results, and (2) to examine cases which we have not analysed in section]^ 
The approach is to determine forces and pressures in the models, using as starting point 
the expressions for the partition function and generating function determined in section 
1^ In the first instance the entropic forces exerted on two confining vertical walls by 
an adsorbing partially directed walk with endpoints tethered in the walls are determined 
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(a) (b) 




Figure 2: The entropic forces exerted by a partially directed walk confined between two vertical walls in 
(as shown in figure]^, (a) The force for walks of length n = 20. endpoint at height h = 0, and with 
a = 1, plotted as a function of w, the distance between the two walls (or equivalently, the width of the 
walk). For small values of w the force is positive (and therefore repulsive, pushing the walls away from one 
another). For large values of w the force is negative; this indicates an attractive force between the two 
walls, (b) Force curves for walks of lengths from n = 20 to n = 200, for a = 1 and h = 0, as a function 
of w. The horizontal length scale is normalised by n for each curve (so that w = [anj for 0 < o: < 1); 
this collapses the curves into a single force curve shown which tends to a limiting curve as n —> oo. Notice 
that the forces vanish at approximately w = |n. 


numerically (see figure [^. By rescaling the forces by the length of the walk, the data 
collapse to a single force curve which is repulsive when the confining walls are close 
together, and attractive when the confining walls are far apart. The typical situation is 
seen in figure when a = 1 and the height of the endpoint of the walk \s h = 0. The 
force curves are modified when a > 1, generally becoming more repulsive as a increases 
(and taking the walk through its transition into its adsorbed phase). 

Secondly, the pressure of the adsorbing partially directed walks on the adsorbing wall 
is examined. If the walk is tethered at the origin, then there is a large entropic pressure 
on the adsorbing wall close to the origin. This pressure decays quickly with distance 
from the origin. If both endpoints of the walk are confined to the adsorbing wall, then 
for walks of length n a secondary pressure peak is seen at a distance about from the 
origin - this peak is the result of the other endpoint of the walk (which exerts pressure 
in the vicinity of the point where the path returns to the adsorbing wall; see for example 
figure 15). 

In section [^we conclude the paper with a few final remarks. 


2. Partially directed walks of width w 

In this section the kernel method (see for example reference ||^|) is used to determine 
G{p) = G{x, y, p, a), the generating function of partially directed walks from the origin 
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Figure 3: The generating function G{p) of partially directed walks from the origin in L^, with edge-visits 
weighted by a and height of endpoint weighted by p, may be classified as above. This classification gives 
a functional recurrence for G{p). Each walk is either (/4) a string of vertical edges (possibly of length 
0), or (6) ends in a horizontal step at height h > 0, or (C) ends in a horizontal step followed by a non¬ 
empty sequence of steps in the up direction, or (D—E) ends in a horizontal step followed by a non-empty 
sequence steps in the South direction (downwards) which does not not step below dh\, or [F — H) ends 
in a horizontal step at height h = 0 (substracted out by H and added back in with weight a in E), or 
{J— K), ends in a horizontal edge at height h = 0, followed by a non-empty sequence of vertical edges, 
subtracted out in K, and then added back with weight a in J. 


in with a collection of activities {x, y, p,, a}, conjugate to length, height of endpoint, 
and number of edge-visits to 9L^. Introduce the following generating variables: 

• horizontal steps: x; 

• vertical steps: y; 

• edge-visits to the adsorbing boundary: a; 

• height of last vertex: p. 


The variable p is the designated catalytic variable, and the generating function of this 
model is denoted by G{p) (the variables {x, y, a} compose a set of parameters or weights 
in the model). The generating function G{p) will be determined by first finding a 
functional recurrence for it, using the classification of partially directed walks illustrated 
in figure 

The contributions to the generating function in figure are labeled by A through 
J, and some are subtracted out while most are added in. Accounting for horizontal and 
vertical steps, and for edge-visits, and height of the last vertex, the following generating 
functions are obtained for each of the terms: 


A 

C 

E 

H 

K 


1 


l-yp’ 

xyp 

1-yp 

xy/p 

1-y/P 


G{p): 

G(y); 


xG(0); 

l^G(O). 


F = ax G(0); 

J = eSG(0); 


Since G{p) = A+B + G + D — E + F — H + J—K, this gives the following functional 
recurrence for G{p)'. 


G{p) = 


1-yp + + l-yp p-y^ 
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Multiply through by (l—yp,)(y — p,) and simplify. The coefficient of G(p,) is the kernel, 
given by 

y(l-(l-x+y2+xy2)^+/i2^ 


K{p) = 


( 10 ) 


(l-y/i)(y-M) 

The numerator of K{p) is a quadratic with roots po and pi such that pi = 1/po and 

( 11 ) 


= ^ (l - X + (1 + x)y2 + ^(1 - x)2 - 2(1 + x2)y2 + (1 + x)^) • 

^ (1 -MMo)(l -MMi) 


Mo 

Notice that 


( 12 ) 


(1 -Aiy)(l -Ai/y) 

SO that AC(0) = 1. 

Series expansion of po and pi shows that pi counts a certain set of walks. Thus, 
Pi is a “physical root" and it will play a key role in determine G{p). The substitution 
X = u{l — v) and y = ^ in F = ^ (pi—y) — 1 gives 


F = 


2uv 


^1 - u{l+v) - \/l- 2u{l + v) + [ 72 ( 1 - v)2j . 


This may be expanded in u to obtain 
F = u + {l + v)u^ + {l+3v+v^)u^ + {l + 6v + 6v^ + v^)u'^ 

+ {l + 10v+20v^ + 10v^ + v'^)u^ H-. 

This is the series expansion of the Narayana generating function given by 


«(-)-eeK:)C-i 

S=1 / = ! V / V 


v'-^uF 


(13) 


(14) 


(15) 


That is, F = N{u, v) and it follows that 

= xy(/\/(x,y) + l)+y. (16) 

The Narayana generating function has useful properties, which will be used later. In 
particular, powers of N{u, v) have nice expansions: 

OO N / \ / \ 

, r / \ / <: \ 

(17) 


(«(o, v)y = J2Yl 

S =1 1=1 

Taking p = 0 \n equation ([^ produces 


s\r+i — lj \i — l 


v'-^uF 


G(y) = 2((l+x(l-^))G(0)-l). 


(18) 


Using equation ( [121 ) and the above in the recurrence (j^ gives 

t 4F G(o) +yp (i+(x.-x-i)G(o))+ iTp - 


i-My ' p-y 
Substituting p = pi gives G(0): 

G(0) = 


(i-y")/^i 


(x(l-y2)(l-a) -y2)/ii +y' 


( 20 ) 
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This is the generating function of walks ending in the adsorbing line (or with endpoint 
at height zero). Substituting this into equation (19) gives a solution for G{p)'. 


(l-MMi)((x(l-y2)(l-a)-y2)/ii +y)' 
Recall that K{pi) = 0; this is useful in simplifying G{p) to 


( 21 ) 


G(m) 


_y-Mi_ 

x{l-ppi){{a-l)pi -ya)' 


( 22 ) 


Substituting x —> xt, and y —> yt, and then doing a series expansion in t, gives the 
series 


1 + {xa+yp)t + {y^ +xy ap+x^ a^ + xy p)F 

+{y^p^+xy^ap^+x^ya^p+x^a^ + 2xy^p^+x^yap+x^yp+xy^)t^ + • • •, (23) 


of walks with t the length generating variable. Substituting pi using the Narayana 
generating function above gives G{p) in terms of A/(x,y): 


X ^_ N[x,y) + 1 _ 

(xa(A/(x, y) + l) - xA/(x, y) - x - l)(xy/j,(A/(x, y) + l) + y/j, - 1)' 


(24) 


Introduce A = x+xA/(x,y) in equation ( [l6| ). Then it follows that pi = y+yA, and 
A = ^pi — 1. Define cr = a—1, then it follows that 


x(1-(jA)( 1 - y/j,(A + l)) ■ 


(25) 


By expanding this in a power series in y, the generating function of walks ending in a 
vertex at height h is 


G/,(a,y) 


Ay"(l + A)" 
x(l —crA) 


(26) 


Substitute A = x+xA/(x,y) and suppress the arguments of N{x,y) = A/ to simplify 
expressions. This gives 


^ , {l + N)y\xN+x+lY 

G,(a.y)= 


(27) 


since a = a —1, and in terms of the Narayana generating function N[x,y). This 
completes the determination of the generating function G/,(a,y) of partially directed 
paths ending in a vertex at height h. 


2.1. The generating function Gy^^{a,y) of partially directed walks of width w 


The generating function of walks of width w and last vertex at height h can be extracted 
from Gh{a,y) in equation ([27|). The strategy is to expand G/,(a,y) in powers of 


A = x+xN, and then to use equation (17) in order to simplify the expressions. 
Expanding the denominator of G/,(a,y) in equation ([2^ in A gives 


G/,(a,y) = i5^c7'^A'^+V'’(l + A)^ 

k=0 


( 28 ) 
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Substituting A = x{l + N) and simplifying leaves 


OO 


Gh{a,y) = ^a^xy^l + Ny 

k=0 


(29) 


Let S{k) be the summand of Gh{p,) above. By expanding the factors in S{k), it follows 
that 

h j+k+1 


S{k) = cj'^xV}^ 


j=0 £=0 






i 


(30) 


so Gh{a,y) = Reverse the order of the summations and allow I to run to 

cx) to simplify the expression. This gives 


OO h 


S{k) = a^x^y^^^i U 

0—r\ i—n / V 


=0 J=0 


h\ fk+J+1 

I 


x^N^. 


(31) 


This introduces powers of Narayana numbers which is replaced by using equation ([Tt]). 
Simplifying the resulting expression gives the following for S(/c): 


OO h 


S{k)=y\^a^}^}_^ 

£=0 J=0 


X 


k+j + 1 


I 


EE 


I 


X’ 

s 

\Fl-\ 


s 

/-I 


X 


1-y' 


y 


2/-2 


(32) 


.5=1 i=l 

The powers of x in the summand is w = k+J+s. This is also the number of horizontal 
steps, and so is the horizontal width w of the walk. Extract the coefficient of x'^ in the 
above to find the generating function of walks of width w and with final vertex at height 
h. This is, after some simplification. 


s{k,w) = J2J2Yl 


I f s 


s 

i+l 


h 

w — k — s 


w-n + l\ y'’+2' 

e ) (T^- 


(33) 


^=0 5=1 /=0 

Here, the summation over s is allowed to run to infinity to simplify the expressions. This 
does not introduce new terms since the binomial coefficients introduce a natural cut-off 
on the sums. It remains to multiplify the above by cr'^ and to sum over k as well - this 
gives the generating function G^hi^, y) of walks of width w ending in a vertex of height 
h, and with a generating steps at height zero. 

A minor simplification can be achieved by noting that the sum over I can be done. 
This reduces the number of summations above to two with the result that 


OO OO , - 

, 1 w + l —s 

s('<.«') = EE—— 


w 


w — k — s) \s — i — \) {1—y'^y 


y 


h+2i 


(34) 


5=1 / = 0 

The 5 = 0 term is equal to y^(„^J, and this should be inserted explicitly in the above. 
Multiplying this by ^ and summing over k gives the generating function for walks of 
width w ending a vertex at height h and with a generating steps at height zero. After 
some simplification, this is 




w+l — s 
w+i+l—s\w—k—s 


s=0 k=0 1=0 

where the subsitution cr = a—1 was made. 


n-1 

i 


w 
s— i 


y 


(7+2// 


a-1)^ 


(l-y2)^ 


( 35 ) 
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The partition function of adsorbing partially directed walks of width w can be extracted 
from Gwh{3,y) in equation (35) by determining the coefficients of (note that the 


coefficient of y^ in Gw/,(a,y) is the partition function of walks of length w+N). In 
other words, it remains only to expand (1—y^)“^ and to collect the coefficient of y'^'. 
Using the binomial theorem, putting n = w + N, substituting cr = a—1 and simplifying, 
gives the partition function of partially directed walks with w horizontal steps and length 


n\ 


n w w 


Zn{w,h) = Y,Y.Yl 


X 


/=0 k=0 s=0 

s —1\ f w \ f h 


w + l — s 
w + i + l — s 


—s 




/ / \s — ij \w — k — sj \^{w + h —n) — s+i 
This expression can be simplified using standard binomial identities to 


n w w 


z,{w,h) = j2J2J2 


w + l —s 


X 


w + / +1 — s 

,=0 k=0 s=0 

s —1\ f w \ f h 


^{n—w — h) — i+s — l 


i J \s — ij \w — k — sj \ s—1 

For example, if n = 8, w = A and h = 2, then 

Zs(4, 2) = 30 + 20 a + 11 3^ + 4 


(a-l)*. 


(36) 


(37) 

(38) 


the partition function of partially directed walks of length 8 and width 4, and with 
generating variable a generating edge-visits at height h = 0. Thus, Zn{w,h) is the 
partition function of adsorbing partially directed walks of length n, width w and endpoint 
at height h (see figure [T]). 


3. The limiting entropic forces 


In this section the entropic forces arising in the limit a —oo are examined. Forces for 
finite values of n will be examined in section (for the non-interacting partially directed 
walk when a = 1). 

The extrinsic free energy of a walk of length n, width w, and endpoint at height h 
is given by 


F„{w. h) = logZ„(w, h) 


(39) 


in terms of the partition function in equation (37). This is the free energy of a walk 
as illustrated in figure [T] confined by two vertical walls a distance w apart. If the walls 
are close together (that is, when w is small), then the walk is constrained, and it loses 
entropy (this reduces its extrinsic free energy). Similarly, if w is large, then the walk is 
stretched in the horizontal direction, and it loses entropy as well. The loss in entropy 
in both these cases induces a restoring entropic force on the opposing walls in figure 
moving them further apart (the repulsive regime), or closer together (the attractive 
regime). At an intermediate distance the entropic force is zero, and this is the zero force 
point. 
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(a) (b) 




Figure 4: (a) The limiting free energy of a partially directed walk as a function of a. The bottom curve 
is the limiting free energy for a = 1 and h = 0. The top curve is the limiting free energy for a = ac and 
h = 0. The derivatives of these curves are the limiting entropic force on the vertical walls, and are plotted 
as a function of a in (b). The darker curve is the force curve for a = 1 and the lighter curve (lower on 
the left and higher on the right in (b)) is the force curve for a = ac. Both curves pass through the zero 
force point at a = |. The curve for a = 1 is symmetric on reflection through the point (|, 0). 


3.1. Entropic forces on confining walls 


The first priority is to determine the limiting (intrinsic) free energy of the walk in the 
limit n —oo, after rescaling edges and lengths in the model by ^ (so that edges have 
length and the vertical walls are a distance w = apart). In this limit the free 

energy can be calculated from equation (35), or more appropriately, by considering the 
exponentially fastest growing terms in equation (37). 

For general values of a and h it follows from equation (37) that 

w+l-s V ^ \/|(n-i/i/-/?)-/+s-l 

W+/ + l-sl / j I 5_/ j I j I 5_]_ 


a — iy < Zn{w, h), 


for any n > 1, for any a > 1, and for any {/, k, s}. Similarly, if k^, Sm} are the values 
of {/, k, s} which maximizes the summand on the right hand side of equation (37), then 

w \ f h \ fk{n-w-h)-im+Sn,-l\ c 


(n+l)^(w+l-Sm) 

W + /'m + 1 — Sm 


1 


W-km-Srr, 


1 


(a-1)^ 


is an upper bound on Zn[w, h). 

Rescale (w, h) by putting w = [ctnj, and h = [cpn]. The summation indices in the 
lower bound are similarly scaled by s = [en\, i = [(5nJ, and k = [Kn\. Taking the power ^ 
and then n —f oo, shows that 


sup 

€,S,K 


2^5^(6 —(5)2(«-'5)(l —Q; —0— 


(40) 


is a lower bound on \\\r\n^oo Zn^'^{[ocn\,[(l)n\). In fact, by the above this limit cannot 
exceed this lower bound for any values of (ct, 0) G [0,1]^. This supremum is particularly 
useful in determining the free energy for the cases that 0 > 0 or a > a^. Determining 
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the free energy for /? = 0 and a = 1, or for /i = 0 and a = ac = 1 + ;^. requires a similar 
approach based on equation(but with h = 0). These cases are considered in sections 
13.1.11 and 13.1.31 below. 


3.1.1. Forces when a = 1 and (p = 0.' Putting a = 1 and /? = 0 in the summand in 


equation (37) reduces the partition function to a single sum. The summand is 

1 ^ + [an\)-[5n\-V 

[5nJ+l 1^(5 


[an\ 
[an\ - [(5nJ 


L(5nJ - 1 

where w = [ctnj and / = [(5nJ. The limiting free energy is obtained by determining the 
maximum exponential rate of growth of this summand in the limit that n ^ oo. This 
can be done by using the Stirling approximation for the binomial coefficients, taking the 
power and then taking a —> oo. This gives 


(a - 2(5)2‘5( 1 + ct - 2(5) 
2“525(a-5)2“(l-cr-2(5)T-“-25)/2- 
The rate of growth of this is a maximum when (5 = (5c and 
5c = \3z ^'>j2oF — 2oL-yl. 

Substituting 5 = 5c gives the limiting free energy as a function of ct: 
Fo(ct) = a log(2ct(l —ct)) — 2ct log('\/C+2ct —1) + log 


(41) 


(42) 


i-Vc 

2a-1 


(43) 


+ (\/C 1) I09 

where C = 2a^ — 2a + l. Notice that Fo(^) = log(l + -\/2); this is not unexpected, since 
partially directed walks grow asymptotically proportionally to (l + -\/2)"+°(”^. 

The free energy is plotted as a function of a in figure Qa). Expanding the free 
energy about its peak at ct = ^ gives 

Fo(a) = log(l + ^/2) - 2V2(a-lf - lV2(a-l)^ + 0((a-l)^), (44) 

where 0((o; —|)®) is a correction which has leading term — ||-\/2(ct —|)®. 

The derivative of JFo(ct) = -^Fo(ct) is the limiting entropic force on the vertical 


a gives the force near a = 


walls confining the walk. Taking the derivative of the approximation in equation (44) to 

1 . 

2' 


= -4V2(a-l) - fV2(a-l)^-^V2(a-l)^ + 0((a-iy). 


(45) 


More generally, taking the derivative of Fo(ct) in equation (43) gives the limiting 


force as a function of a. This is plotted in figure |^b) and simplifies to 


T'o(ct) = log 


2Q!(l-a) 


+ ^y^log 


(C-a)( 2 a+C-lf 

(C+a)(l-CF 


(46) 


(C + 2q:-1)2; ' 2C 

It may be verified that To(cf) = —JFo(l —ct), so JFo(|) = 0. That is, the limiting force 
curve is symmetric on reflection through the point (|,0). For small ct > 0 the above 
shows that T'o(ct) 


logct|; this will also be shown later (see equation (65)). 
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Figure 5: (a) The limiting free energy 4>) plotted as a function of a for d e {0, 0.1, 0.2, 0.3, 0.4, 0.5}. 
Notice that a G [0,1 —(/>]. (b) Force curves tF^ia.d) as a function of a. Increasing d moves the zero 
force point to smaller values of a. For d = 0 the zero force point is located at 0:^(0) = j. 


sup 

€,(5 


3.1.2. Forces when a = 1 and (p > 0; The limiting free energy can be determined by 
computing the supremum 

a“(/)‘^(l-a-(^-2(5+2e)(i-“-‘^-26+26)/2 1 

which is found after putting k = 0 and by taking a —f 1“ in equation (40). The 
supremum is realised when e = Cs, and 5 = 5s, where Cs and 5s are given by (note that 
0 < ct + 0 < 1) 

c^{c^ — c^'^ + ^y{2a^ — 2a + l)(p'^ — (p^) 


5s = 


{a+(p)[l+(p-a) 

a {l — a — (p){(p+(p^ — ^y{2a'^ — 2a + l)(p'^ — (p^) 


2<p 


{a+(p){l + (p-a) 


(48) 


(49) 


Substituting these values to determine the supremum (using Maple ||^|) gives a very 
lengthy expression for the free energy F^{oc,(l)) as a function of (q:, 0). It is plotted 
against ct for various values of (p in figure [^a). 

The force iF(i>{a, (p) is the (partial) derivative of Fcp{oc, (p) to a. This is plotted 
against a for various values of cp in figure [^b). If 0 = 0 then the force vanishes when 
ct = ^, so that T'^(|, 0) = 0 (where cp = 0). This is the zero force point for the case 
0 = 0. The zero force point is generally located at adcp), so a<-(0) = |. If 0 increases, 
then the zero force point moves to smaller values of ct. By expanding the force in 0 and 
a, the location of the zero force point as a function of 0 can be estimated by determining 
the first few terms in a series expansion in 0. To order 0® this is 

c^c(0) = \-\V2(P^ + |(3V2-2) 0^ + ^(144-14972) 0® + 0(0®). (50) 

This estimate is surprisingly accurate. For example, ac(O.l) = 0.4964... and the 
estimate is ac(O.l) 0.4965..., and ac(0.3) = 0.4669... while the estimate is 
ac(0.3) 0.4717.... The estimate for 0 = 0.5 is ac(0.5) 0.4249... as opposed 

to the exact value ac(0.5) = 0.40095.... 
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3.1.3. Forces when a = ac and (p = 0; Put a = ac = 1 + :^ summand in equation 


(37) and bound the partition function in a way similar to the way it was done in section 
(3.l|). Take the power and let n —> oo. This shows that the limiting free energy can 
be determined by computing critical values of (e, <5, k) to find the supremum 

a“(l-a-25+2e)(^-“-25+24/2(^)« 


sup 

€,S,K 




(51) 


The critical values of (e, 5, k) in the above are 


Cs = 4(2 —V2) f l + a +Y 4 'y 2 o;(l —o;) + (l + a)2 


5s =k2-V2)(l-cr); 


Ks = i(2-V2) (l + ^/2)2cr-l-^/4^/2a(l-a) + (l + cr)2 


Substituting these, simplifying, and then taking the logarithm gives a lengthy expression 
for the limiting free energy Fc(cr) as a function of a. Expanding this in a about ct = | 
gives, after simplification, 

Fc{a) = log(l + y2) - 1(272-1) {a-lf + ^(72-1) (o;-|)^ 

r(l + 72)7l-272)7l +872) + 0((cr-|)^). 


2 y5^-L| ^-L ^ V ^ J y J. \ KJ v J 2^ I V, 2 ' (52) 

Notice that Fc(ct) is not symmetric about ct = |, unlike Fo{a) (the free energy when 
0 = 0 and a = 1). Taking the derivative of Fc(cf) to ct gives the force near ct = F 

^,(a) = f ( 2 V 2 -I) (a-i) + 3^(V2-l) (a-l)2 

- 32^(l + V2)2(l-2V2)'(l+8V2)(a-i)=' + 0((a-i)‘'). (53) 

The full expression of the free energy Fc{oc) is a lengthy expression determined using 
Maple and it is plotted in figure|^a). Taking its derivative to a shows that the force 
vanishes when ct = |, consistent with the series expansion above. Thus, JFc(|) = 0. 
The force curve is plotted in figure Qb). 


3.1.4. Forces when a = ac and (p > 0; If a = a^ and 0 > 0, then the critical values of 


(e, 5, ac) in equation (40) are given by 


£5 = 1-2- 


|(2 —72) ^1 + ct —0+^472 q:( 1 —ct —0) + (1 + q: —0)2 Y 


5s = 1(2-72) (l-cr-0); 


ACs = T 2- 


^(2-72) (^{l + V2ya-3(p-l-^4V2a{l-a-(p) + {l+a-(py], 


Substituting these into equation (40), simplifying, and then taking the logarithm, gives 
a lengthy expression for the limiting free energy Fc(ct, 0) as a function of a and 0. The 
entropic forces are again given by the partial derivative JFc(o;, 0) = ^Fc{oc, 0). 

In figure [^the limiting free energy and forces are plotted as a function of ct for fixed 
0. With increasing 0 the zero force point ctc(0) moves to smaller values of a. The 
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(a) (b) 




Figure 6: (a) The limiting free energy Fc(a) plotted as a function of a for 0 G {0, 0.1, 0.2, 0.3, 0.4, 0.5}. 
Notice that a G [0,1 —(/)]. (b) Force curves T'c(a) as a function of a. Increasing (p moves the zero force 
point to smaller values of a, and this point is located at exactly | — ^0. 


location of this point may be approximated by expanding the force in 0 and in ct, and 
then solving for ct as a function of 0. The result is 

OLci(p) = ^ - ^(p + 0{(p^). (54) 

Numerical analysis of the model suggests that the correction 0(0^) is zero, so ctc(0) is 
equal to | —This is verified by substitution and then simplification of the expression 
for the force. That is, the zero force point is located exactly at 

a^c(0) = (55) 

if a = ac. 


3.1.5. Forces when a > ac and (p > 0.' The adsorbed phase of the model is obtained if 
a > ac. In this case the critical values of (e, <5, k) in equation ( |40| ) are given by 

65 =^ ol — (p-\- \/S{a — l)o;(l — o; — 0) + (l + o; — (p)^ 

5s = ^ (l-a-0); 

Ks = :^ (^{^a-l)oL-3(p-l-\/8{a-l)a{l-a-(p) + {l+a-(py 
provided that a > a^. 

Notice that Gs + Ks = a — ^cp < ol as required in equation (j^), provided that act > (p. 
Similarly, cp—oc+es + Ks = 0(1 — ^) > 0, provided that a > 1. In addition, 65 > 5s and 
l — a — (p — 26s > 0 if both a > 1 and a+(p < 1 . 

Substitution of the above critical values for (e, 5, ac) and taking the logarithm gives 
the limiting (intensive) free energy F{oc,(p). This expression for F{oL,(p) is lengthy (it 
was determined using Maple |[^) and will not be reproduced here. 

Taking the derivative of F(ct,0) to ol gives the limiting entropic force lF{oL,(p). 
Putting both 0 = 0 and a = 1 simplifies lF{oL,(p) to equation (46), and this force is 
plotted in figure l^a). 
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(a) (b) 




Figure 7: (a) The entropic force curve for a = 1 and (^ = 0 is recovered from the force T'(a,0). This 
curve is identitical to the force curve tFo for a = 1, plotted in figure Qb). (b) The force curves tF{a, a) 
for increasing values of a. If a = a^, then the zero force point is at a = Increasing a moves the zero 
force point to the right and flattens the force curve, compared to the force curve JTq plotted in (a). 


In figure [^b) the forces for values of 0 = 0 and a > and for 0 = 0 are plotted. 
The expressions for these forces are lengthy, and were obtained using Maple (^|. At the 
critical adsorption point ac = 1 + :^ bhe zero force point is located at ct = |, unchanged 
from its location when a = 1. 

Increasing the value of a gives force curves which are progressively higher in the 
graph (and so larger if repulsive, and less attractive, if attractive), while the zero force 
point moves towards larger values of a. If 0 = 0, then the location of the zero force 
point may be estimated by expanding the force in a and solving for the zero force point. 
This gives the following asymptotic expression for the location of the zero-force point: 
If a —ac is small, then the zero force point is located at 

OLc[a) = I + |(a-ac) - |(\/2-l)(a-ac)^ + 0((a-ac)^). (56) 

For example, if a = ac + 0.1, then ctc(ac + 0.1) = 0.5488 ... while the approximation gives 
OLc[ac + 0.1) se 0.544.... Similarly, ac(ac + 0.05) = 0.5247... while the approximation 
gives ctc(ac + 0.05) ss 0.524 .... 

For large values of a an asymptotic expression for the zero force point may be found. 
This is given by 

= 1 — (3^1)3 + (a-l)4 “ (a-l)5 +0[[a—l) ®). (57) 

This becomes very accurate for large values of a. Numerically ctc(4.0) = 0.9559 . .. while 
the above approximation gives ctc(4.0) ss 0.934 .... Similarly, ctc(6.0) = 0.9884... while 
the approximation gives ctc(6.0) ss 0.987.... 

4. Finite size asymptotics for a = 1 and h = 0 

In this section an asymptotic expression for the partition function in equation ( [37| ) (this 
is for the special case that a = 1 and /? = 0) is determined. Notice that if /? = 0 then 
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only terms with w — k — s = 0 survive, and if a = 1, then k = 0, so that w = s. This 
reduces the triple summation in equation (37) to a single summation: 

w—l\/' w \ /'^{n + w) — i—1 


Z,{w, 0) = 


W 


^ / + ! 

/=0 


/ 


w 

w — i 


w— 1 


(58) 


This partition function will be approximated by an integral, and the integral will be 
approximated by a saddle point formula. The first step is to find a good approximation 
for the summand. 

Thus, consider the summand in equation (^) and denote it by S. Its binomial 
coefficients will be approximated by using the Stirling approximation for the factorial: 


n! = V2'Knn'^e 


'(1 “ 1^ + 2M7T + (59) 

Take the logarithm of the summand in equation ( |4^ and expand the binomial coefficients 
into factorials. This gives 

log S = log(/+l) — log(M/ —/) + log((|(n +w) —/ —1)!) 

- log((|(n + i/i/) —/ —1)!) — 2 log((M/ —/ —1)!) + log(i/i/!) — 2 log((/+l)!). (60) 

Substitute the factorial by the Stirling approximation above, and simplify. This gives an 
expression approximating the logarithm of the summand. 

In order to implement a saddle point approximation, it is necessary to find the 
location of the saddlepoint in the summand. Thus, in the simplified approximation of 
log S above, substitute n = ^, I = ^ and w = |. Expand the resulting expression in e 
and collect the leading term (which is 0(e“^)). This is a lengthy expression in 5 and a, 
and the saddle point is located by taking its derivative with respect to 5 and solving for 
the stationary point. This gives 5 = | ± \\/2oF — 2olF\ as the possible locations of the 
saddle point(s) (see for example equation (|4^). Notice that the choice of the plus sign 
before the square root gives 5 > ct - so that / > w in the summand. This is outside the 
range of / in equation (44), with the result that the minus sign before the square root 
is the correct choice of the sign. Thus, the saddle point is located at asymptotic values 
of / in the summand where / = [5sn\, where 

5s = 5 - 5V2a2-2a + l. (61) 

Numerical work shows that the width in the peak about 5s is proportional to y/n, 
thus the summation in equation (44) will be approximated by an integral by putting 
^ % and / = ^ + f in the logarithm of the summand. The result is expanded 

in a Laurent series in e; this gives leading terms of order O(^) and 0(\). The series 
is truncated to 0(e), then exponentiated and simplified, before it is integrated over 5. 
The result will be a saddle point approximation to the partition function Zn{w,0) with 
w = [ctnj, once the substitution e = ;^ is made. After significant symbolic computations 
calculations using Maple (^, the approximation to Zn(M/, 0) for a = 1 is 

Z„(LanJ ,0)-(62) 


7rn2cf(ct —l)2y {l — a)^/C + C 

1 - 3ct + 3^2 + {l-2a)^/C' 
o;“(l-o;)“ 


'yc + 


a 


nl2 


\fc — OL 


X 














Forces and pressures in adsorbing partially directed walks 


18 



Figure 8: The exact and approximate free energies Fn{w,0) for n = 20 (bottom curves) and n = 100 
(top curves). The approximate expression is obtained by Fn{w,0) = \ogZn{w,0) where Zn(M/,0) is 


approximated by equation (62). On this scale there is little differences between the exact curve and its 
approximation. 


where C = 2oF — 2ol-F1. 

The approximation for above also approximates the extensive free energy, given 
by Fn([ctr!j,0) = log Zn([ctnJ, 0) according to equation ( [3^ . In figure the exact 


extensive free energy, and its approximation (by the logarithm of the right hand side 


of equation (62)), are plotted on the same graph as a function of - for n = 20 and 


n = 100. On the scale of this plot, the approximation is indistinguishable from the exact 
curve. 

An asymptotic formula for the entropic force is obtained by taking the derivative 
of log Zn([aJ, 0) in equation (62) to a while recalling that w = [ctnj. In figure the 
exact force curve, and the asymptotic approximation to it, are plotted on the same 
graph for n = 200. This shows that the forces are very small for ol = \. In fact. 


taking the logarithm and then the derivative of the right hand side of equation (62), 


gives a complicated expression which simplifies to 1 + V2 if ct = |. That is, there is no 
dependence on n at this point, and it may be verified, using symbolic computations (^ , 
that the coefficient of n vanishes at this point. That is, the asymptotic expression of 
the finite size force is of the form 


^ Aa + B^n 


and Aa and Ba are functions of a such that A 


1/2 


(63) 

= 1 + v^ and Bi /2 = 0. Thus the force 


does not vanish at ct = ^, but is, instead, a constant value, independent of n, where 


T'n(l) = I + V 2 . Since ^'^(o;) increases with n for ct 7 ^ this indicates that 

n ^ ^ 1 

Zn —> 0 , as n ^ 00 and a ^ F 
Tn{a) ^ 2 


(64) 


One may similarly use equation (62) to determine an asymptotic expression for the 


entropic force on the two vertical walls in figure [T] This may be done by taking the 
derivative of the free energy to ct. For ol small the result is 

Tn{oL) ~ + {n-\) - nioga + 0{na). 




(65) 
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Figure 9: The normalised force curve for n = 200 and a = 1 (see figure |^b)) plotted as a 

function of w. There are two curves, namely, the exact force for n = 200, computed from equation (66), 
and the asymptotic approximation to tFnict) for n = 200. The curves pass very close to | when tFn = 0. 


Fixing w = [nctj and then taking ct —f 0+ shows that 
values of n. Moreover, ~ | logctl (see section 3.1.1). 


n\ log ct| —f oo for small 


5. Numerical calculation of the entropic forces 

5.1. Entropic force on the vertical wall 

For finite values of n and w the entropic force may be defined by the finite difference 

h) = ^ {Fn{w+2, h) - Fn{w, h)) (66) 

of the extensive free energy. The difference between the free energies for widths w+2 
and w are used in order to avoid parity effects in the model; these parity effects are due 
to fact that the underlying square lattice is a bipartite graph. 


5.1.1. The case a = 1: In this case there is no attraction between the walk and the 
adsorbing line, and the model is in its desorbed phase. 

In figure a) the entropic force T' 2 o('^. 0) is interpolated on discrete points as a 
function of the width w for walks of length n = 20 and endpoint fixed at height h = 0. 
As expected, the force between the confining walls is large positive (repulsive) for small 
values of w, and large negative (attractive) for large values of w. At a critical value of 
w, namely w 'ss 10, the force is zero, neither attractive nor repulsive. 

Forces of paths of different lengths n can be compared by scaling w by n in figure 
[^a). In figurej^b) the interpolated force curves for n G {20, 40, 60, 80,100} are plotted 
against a = ^. These curves collapse to a single curve, namely a universal force curve 
which is the limiting force in the limit n —> oo with length rescaled by n. 

Note that the force curves in figure]^ pass very close to the point {w, n) = {^n, n) 
when lFn{w, 0) = 0. That is, the force is approximately zero when w = ^n, or when the 
confining walls are w = ^n apart. This occurs when the walk gives ^n horizontal steps, 
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Figure 10: The entropic forces exerted by a partially directed walk between two vertical walls. The width 
w of the walk (or distance between the confining walls) is rescaled by n (the length of the walk), (a) The 
entropic forces for n between 20 and 200 plotted for walks with a = 1 and height of endpoint h = \n. 
Notice that this distorts the curves when ^ is larger than about 0.5. (b) Similar to (a), but now with 
h = |n. 


and this is consistent with the location of the zero force point at ct = | determined in 
section 13.1.11 

The convergence of the finite size forces to the limiting force curve is quick, this 
can be seen, for example, in figure where the limiting force, and the exact force for 
n = 200, are plotted on the same graph. These curves are indistinguishable on this scale. 

In figure [T^ the effects of the height of the endpoint on the induced forces are 
displayed. In figure 10(a) the force curves are plotted for h = [cpn] with 0 = |- Since w 
is rescaled by n, the curves again collapse to an underlying force curve. For small values 
of w the force curves are not changed significantly from those in figure [^b). In figure 
[l^b) similar data are presented, but now with h = |n. 

Finally, the effects of the adsorption activity a is examined in figures [II|(a) and 
0 :b)- For a larger than 1 the walks are attracted to the adsorbing boundary 9L^, and 


they adsorb at a = ac = 1 + ^- In figure 11(a) the forces are plotted for a = 2 > a 


The data shows a longer range of repulsiveTorces, and for large w, when the forces are 
attractive, the magnitude of the forces are reduced. This distortion of the forces can 
be seen even in the case that a = ac, as shown in figure [T^b), although here it is more 
subtle. There are, however, an apparent reduction in the magnitude of the attractive 
forces when w approaches n. 


5.2. The pressure on the adsorbing wall 

Let u = ((w, 0 )~(m/ + 1, 0)) be an edge in the adsorbing boundary a distance w 
from the origin. The pressure Fln(i/i/) on u is the change in extensive free energy if u is 
a blocked so that walks cannot pass through it. That is, if T is the free energy of walks 
of length n, and fn(a) is the free energy of walks avoiding the edge u, then the pressure 
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(a) 


(b) 



Figure 11: The entropic forces exerted by a partially directed walk between two vertical walls. The width 
w of the walk (or distance between the confining walls) is rescaled by n (the length of the walk) and the 
height of the endpoint of the walk is at h = 0. (a) The entropic forces plotted for a = 2. The force curves 
are shifted in the positive direction compared to the results in figurej^b). This shows that the forces are 
repulsive for larger values of w, and, when they become attractive (negative), then they are weakened, 
(b) The entropic forces plotted for a = ac (the critical adsorption point in the model). 


on u is given by 

Un{w) = fn- fn{u). (67) 

The (extensive) free energies F and fn{u) are determined by the corresponding partition 
functions, and the partition function of walks avoiding u will be found by first computing 
the partition function of walks passing through u. 

As in equation (^), denote the partition function of partially directed walks of 
length n from the origin in L|, with edge-visits to weighted by a, of horizontal 
width w and with endpoint at height h, by Zn{w, h). This partition function can be used 
to determine the partition function of walks passing through an edge u as illustrated in 
figure 

Let bn{v, h) be the number of partially directed walks of length n in L^, from the 
origin, making v edge-visits to 9L^, and with final vertex at vertical height h. The 
partition function of this walk is given by 

n 

Bn{h) = '^b,{v, h)a\ (68) 

v=0 

and it is implicitly a function of the adsorbing activity a. If p is introduced as the 
generating variable conjugate to h, and t is conjugate to n, then the generating function 
of this walk is given by 

t[l — pt) (l — [l+2p)t + F + t^ — v/D) 


6 (t, p, a) = 


(69) 


((l + p+p'^)t - {l + F + t^)p) ((l + a)t2 + (i_a)(l-t-Ft3-A/D)) 

Summing Bn{h) over h gives the partition function of walks ending at any height, namely 
Bn — Yl!h=o Bn{h), and this is also implicitly a function of a. 
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1 / 1 / = [an\ 


Figure 12: A partially directed walk passing through an edge in the adsorbing boundary a distance \n 
from the origin may be calculated by concatenating a bargraph path of width w with a partially directed 
walk as shown. The bargraph path and partially directed walk are joined by inserting an edge-visit 
u = ((w, 0)~(w-F 1, 0)) between the endpoint of the bargraph path and the first vertex in the partially 
directed walk. Since u is an edge-visit, it carries weight a, similar to other edge-visits in the model. 


Notice that in terms of Zn{w, h), 

n n n—w 

B„(h) = Yi Z„(w. h). and P« = Y'Z, 


(70) 


w=0 w=0h=0 

The partition function of adsorbing bargraph paths is 6^(0). 

The partition function of walks passing through the edge u may be obtained 
concatenating a bargraph path of width w, with a partially directed walk ending in a 
vertex at height h. That is, concatenate a walk generated by Z^(m/, 0) with a path 


generated by as illustrated in figure 12 and then sum over i. That is, the 


partition function of walks passing through the edge u and ending in a vertex at height 
h is 


n-l 


= aJ2z^{w,0) Bn-e-i{h). 


(71) 


l=w 


Similarly, for walks ending in a vertex at any height, the partition function is given by 


n-l 


vy = aYZtC^.O) 


(72) 


e=w 


Note that the extra factor a in these expressions is the weight of the edge u, which is not 


accounted for in the partition function otherwise. See figurefor more explanation. 

The pressure on the edge u by walks of length n can be found from equation 
(67) by determining fn and fn{u). In particular, for walks ending a vertex at height h, 
fn = log Bn{h), and fn{u) = \og{Bn{h) — Yn^\h)) since the partition function of walks 
ending in a vertex at height h which avoids the edge u is Bn{h) — Yn^\h). Thus the 
pressure on u by walks ending in a vertex at height h is given by 


n!.'>( 


W 


, b) = log Br,{b) - log {Bn{h)-Yy\h)) = - log (^1 


Bnih) 


(73) 
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Figure 13: The pressures and as a function of w. Notice that w is normalised by 

n and that the pressures are scaled by multiplication with n. In these simulations a = 1 and n = 100. 
The pressure of a bargraph path n(^^(u/,0) on an edge u in the adsorbing line at first decreases quickly 
with increasing w, but then goes through a secondary, peak strongly modified by a parity effect, at about 
w = \n. For w approaching n the pressure decreases to zero. The pressure of partially directed walks 
with unrestricted endpoint decreases monotonically with increasing w to small values for w > \r\. 


Next, the 


pressure 
'(!)/ 


0) in equation (73) can be approximated by 
approximations of V)}^•'(0) (see equation (71)) and the bargraph partition function 6n(0). 

is approximated by using the asymptotic approximation (62) for Zn(w,0), and 
for 6n(0) the asymptotic approximation 


Bn{0) 


y/Trrd 


r-\ '1+3/2 

l + ^/2) 


1 - 


36+21^/2 , 1745+1260\/2 , ^,„-3 

I occ „2 “T w{n 


16n 


256 n2 


))■ 


(74) 


derived in reference |[T^, can be used. Substitute these in equation (|7p, and approximate 
the summation by an integral (notice that the summation in equation ( [tI] ) includes values 
of i close to zero, and again close to n, where the approximations in equations (62) and 
( [74l ) are poor). Thus, by restricting the summation to w > A, the resulting approximation 
will be 

n—5 

yW( 0 )~a z,(w, 0 )e„_,_i( 0 ). ( 75 ) 

i=w+l 

That is, using this approximation, and the approximation for 6^(0) above, nn^^(M/, 0) 
may be approximated using equation (73). This approximation is interpolated on odd 
values of w for n = 100 in figure (there are strong parity effects, as seen in figure 


13 ). 


Similarly, the pressure of walk ending in a vertex at any height on the edge u is given 


by 


= log P„ - log (P„-yy) = - log (l-T"\ 


(76) 
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Figure 14: The pressure of bargraph paths of length n = 100 and a = 1 on dh\ interpolated 

on odd values of w. The pressure is large close to the origin, declines steeply with increasing w, and then 
passes through a secondary peak before decaying to zero. The estimated pressure were computed by using 
the approximation in equations (74) and (75) in equation (73). 


since fn = Pn in equation (67) (where Pn is the partition function of partially directed 
walks of length n ending at any height), and fn{u) = Pn — yn^^ is the partition function 
of walks of length n ending at any height, and avoiding the edge u. 

The pressures h) and r\^i^\w) can be computed for small values of n (say 

n < 100) and involves quadruple summations (three in equation (37) and another in 
equation (70)). These expressions were coded into a C program in order to explore 
the pressures and Numerical simulations for a = 1 show a strong parity 

dependence on w (the distance of the edge u from the origin). For example, in figure 


13 the pressures are displayed as a function of w (normalised by n) for n = 100. The 
pressure 0) decreases quickly with increasing w, but shows a secondary peak for 

w approaching \n, with a strong correction due to parity effects. The pressure r\^n\w), 
for walks with endpoint at any height, decreases monotonically to zero with increasing 
w. Parity effects for small values of n remains visible in the results; thus, in order to 
suppress these, pressures were evaluated only for odd values of w, in what follows below. 

In figure ^(a) the pressure Fln^^w, 0) for bargraph walks are interpolated on odd 
values of w and plotted against for n G {20,40,60,80,100}. For any fixed n the 
pressure decreases with increasing w, but there is a secondary peak in the pressure when 
w approaches ^n. The appearance of this secondary peak may be understood by noting 
that the endpoint of bargraph paths in would cluster around the expected horizontal 
width of the path, and so exerts pressure on edges in this vicinity. 

It is also noticeable that for fixed values of ^ the pressure declines with increasing 
n. This follows because, with increasing n, the likelihood of a walk passing through 
the edge u = (m/~(m/ + 1)) decreases. Since the rescaling of w by n effectively gives 
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(a) (b) 




Figure 15: The pressures as a function of w with a = 1. (a) The pressures for n in 

{20,40,60,80,100} against With increasing n the pressure decreases, but it shows a secondary 
peak when w approaches |n. (b) Since the length of a rescaled edge decreases proportional to ^ with 
increasing n, the pressures in (a) can be rescaled by plotting nfln^^ against This brings the data in (a) 
closer together and suggests an underlying limiting rescaled pressure. Note that the data in these graphs 
are interpolated on odd values of w. 


(a) 




Figure 16: The pressures as a function of w with a = 1. (a) The pressures for n in 

{20,40,60,80,100} as interpolated as a function of With increasing n the pressure decreases 

monotonically. (b) Since the length of a rescaled edge decreases proportional to ^ with increasing n, 

the pressures in (a) can be rescaled by plotting nFIn^^ against This collapses the data in (a) closer 

together and suggests an underlying limiting rescaled pressure. Note that the data in these graphs are 
interpolated on odd values of w. 


edges of length a multiplication of 0) by n should rescale the pressure curves 

to compensate; there may even be a limiting (non-zero) pressure curve in the limit as 
n —> oo in this model. This possibility is explored in figure[I^b), and while there remains 


some spread in t 
In figure 16 


le curves, they do cluster nicely together. 

a) the data for the pressure is interpolated on Since the 

endpoint of the walk is unconstrained in this case, the pressure declines to zero with 
increasing w without the presence of a secondary peak (see figure [T^a)). These data 
can also be rescaled by plotting in figure 16'b) against -. This exposes, as in 
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(a) (b) 




Figure 17: The pressures 0) (for partially directed walks with endpoint at height h = 0 and with 

3 = 2) and (for partially directed walks with endpoint at any height h and with a = 2) as a 

function of w. (a) The pressures forpartially directed paths walks (or bargraph paths) with endpoint at 
height h = 0 for n in {20,40, 60, 80,100} as a interpolated as a function of For w less than about 
the pressures are large, and while slowly decreasing, may approach a non-zero constant as n —> oo. For 
w large the pressure approaches zero. Note that there is no-rescaling of n(^)(n/,0) by n in this figure, 
unlike the graph displayed in figure [T^b). (b) The pressures for walks with endpoint at any height, for 
n in (20,40,60,80,100), as interpolated as a function of These pressures were also not rescaled by 
n, unlike the data in figure [T^b). For small values of w the pressures are large and non-vanishing, but 
for large values of w the pressures approaches zero. Note that the data in this graphs are interpolated on 
odd values of w. 


onto dh\. 


figure [T^b), an apparent limiting pressure curve. 

n figure[^the pressures are displayed for a = 2; that is, when the walk is adsorbed 
Comparison to the data in figures 15'a) and 16[a) shows that the pressures 
here are larger, and moreover, for small values of w, do not decrease markedly with 
increasing n. In fact, for values of w less than about \n both the data for and fln^^ 
suggest that the pressures converges to a non-zero value as n becomes very large. For 
large values of w, the pressures approach zero. 


6. Conclusions 


In this paper a two dimensional partially directed walk model was used to examine the 
forces and pressures of a linear polymer near an adsorbing wall. We determined the 
generating function of the model, and used it to compute the entropic forces it induces 
between two vertical confining walls, and also to determine the pressure exerted on the 
adsorbing wall. 

Similar lattice models of the entropic forces in polymers have been examined 
elsewhere. For example, in reference |[^| a self-interacting self-avoiding walk model 
of pulled walks was examined in LF using exact enumeration data. In the constant 
distance ensemble the force-extension curves were determined for walks below, at, and 


above the 0-temperature. The results (see for example figure 9 in |13|) show that the 
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force is attractive (pulling the end-points of the walks together) if the distance between 
the endpoints is a non-zero constant fraction of n (the length of the walk). This is not 
surprising, since the natural length scale in this model is n‘^ (where u < 1 \s the metric 
exponent of the self-avoiding walk). Thus, there is no zero force point in this model. 

It is possible to prove that there is no zero force point in a model of self-avoiding 
walks confined to a slit of fixed width w. This may be done as follows. A bridge is a 
self-avoiding walk from the origin, with first step in the x-direction, never to return to 
the X = 0 plane, and with last vertex having maximal x-coordinate (see, for example. 


reference |[^ for a definition). If a bridge has x-span equal to w, then it has width w. 
If the number of bridges of length n and width w is denoted /3n(i/i/), then the growth 
constant is known to exist (see, for example, reference |[^|), and is given by 

\ogp^= lim ^\ogpn{w). (77) 


The force on the vertical walls is given by the discrete derivative 
Fw log log 


(78) 


It is known that Pw-i < Fw < where p is the growth constant of the self-avoiding 
walk |[T^|. This shows that > 0 for all values of w, and there is no zero force point 
in this model. Moreover, taking w —oo shows that lim^z-^oo F^ = 0, since p^ ^ p as 
w ^ oo (T^. 

An alternative model, with length in the horizontal direction scaled by n is obtained 
by choosing w = [ctnj in the above model. The partition function of this model if 
defined by Zn{oL) = /3n([ctnJ). The (extensive) free energy of this model is given by is 
Fn{oc) = log Zn{oL). The microcanonical density conjugate to the width of bridges in this 
model is defined by the limit 

P{a) = lim iPniloLn]))^^" . (79) 


It can be shown that this limit exists (see reference Il5|), and is a log-concave function 


of ct G [0,1], and so is differentiable for almost every a G (0,1). 

Define the free energy of pulled bridges in the constant force ensemble by 




lim 


log '^Pniw) 


X' 


(80) 


w=Q 


where x = e^ (and f is the horizontal pulling force on the endpoint of the bridge). This 
limit can also be shown to exist (see reference |[^| for more details). The microcanonical 
density function P(ct) is related to the Legendre transform of T'(x), namely 

log P(ct) = inf {T’(x) — ct log x}. (81) 

x>0 

It is known that tF{x) is a convex function of logx and that 




(82) 


= log p, if X < 1; 

> log p, if X > 1. 

For proofs of these facts, see references [l[[^. This implies that, (1) the right derivative 
of P(ct) at ct = 0 is equal to 0, and (2) P{oc) is strictly decreasing for a G (0,1]. 
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The limiting force on the vertical walls in this model is given by 

Fa = ^ log P{oc). (83) 

By the properties of P(ct), this shows that is a strictly negative force for all values of 
a e (0,1) (implying that it pulls the vertical walls together). Moreover, Fq = 0 and Fa 
is a (strictly) decreasing function of ct > 0, showing that the magnitude of the pulling 
force increases as the endpoints of the bridge are taken further apart. That is, there are 
no zero force points in this model as well. 

In contrast to the results for the self-avoiding walk discussed above, our results 
show that the forces between confining walls may be either attractive (if the endpoints 
of the walk are far apart), or repulsive (if the endpoints are close together). Generally 
the strength of the repulsive forces increases if the walk adsorbs on the adsorbing wall, 
as seen, for example, in figure [T^b). 

A particular interesting result is the location of the zero force point in the models. 
For bargraph paths (when h = 0 and a = 1 in section |3.1.lD this point is located at 


q:c( 0) = so that the force vanishes when the horizontal extent of the path is one-half 
its total length. The location of this point does not move with increasing attraction into 


the adsorbing line in the desorbed phase, as seen in section |3.1.31 where at the critical 
point ac, the zero force point is still located at ct = ^. In other words, for all desorbed 
bargraph paths, the zero force point is located at ct = | (this may be verified by using 
the full expression for the free energy F(a,0) in section 3.1.5 (and by putting 0 = 0 
therein). 

For values of a > ac and 0 > 0 the location of the zero force point becomes more 
interesting. In particular, for a = ac the location of this point is a function of 0 given in 


equation (55), and approximate locations for the zero force point in the adsorbed phase 


(where a > ac) are given for 0 = 0 in equations ( [56D and ((STj). 

We also succeeded in developing an asymptotic formula for the partition function of 
the partially directed walk in the case that a = 1 and h = 0 (these are bargraph paths). 
This result was used to verify our results, and it shows excellent agreement. 

We examined the pressure of the walk on the adsorbing wall - this is always positive, 
and is sharply peaked near the origin (where the walk is tethered to the adsorbing wall). 


Rescaling the pressure (see for example figure 15) shows a secondary peak at the average 


position of the freely moving endpoint of the walk. 
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